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At metal/insulator interfaces, the metallic wave functions with an energy eigenvalue within the 
band gap decay exponentially inside the insulator, creating a continuum of metal-induced gap states 
(MIGS). Despite being a genuine interface property, the features of MIGS have been typically dis- 
cussed in terms of bulk Bloch functions of the insulator with an associated complex wave vector. 
We have developed a model to compare the energy-dependent decay factors obtained from the first- 
principles complex band structure of a bulk insulator (PbTiOa) and those coming from simulations 
of realistic capacitors, including explicitly the interface with the electrodes (Pt and SrRuOs). The 
agreement is very good, especially for the noble metal-based capacitor. For metal oxide electrodes, 
symmetry filtering arguments are required to understand the fine details. The ferroelectric polar- 
ization within the insulator layer tends to increase the effective decay factor, implying MIGS with 
shorter penetration lengths. These results are important in the design of giant tunneling electrore- 
sistance and magnetoresistance devices with metal/oxide junctions. 

PACS numbers: 73.40.Gk, 73.30.+y, 77.80. Fm, 31.15.A- 



Intcrfaces between different oxide materials exhibit a 
fantastic variety of functional properties, some of them 
intrinsic to the boundary between the constituent com- 
pounds [1] . A ubiquitous type of interface is that formed 
between a metal and an insulator or semiconductor, since 
the presence of electrodes is required for the applica- 
tion of electric fields in solid state devices. At any 
metal/insulator interface, bulk Bloch states in the metal 
side of the junction with eigenvalues below the Fermi en- 
ergy and within the energy gap of the insulator cannot 
propagate into the insulating side. These states do not 
vanish right at the interface cither, but they decay expo- 
nentially as they penetrate into the insulator [2] . Indeed, 
the tails of the metal Bloch states might have a significant 
amplitude for a few layers from the interface, creating a 
continuum of gap states [the so-called metal-induced gap 
states (MIGS)] [2]. These MIGS are essential to deter- 
mine many interfacial properties. The transfer of charge 
associated to them contribute to the interfacial dipolc, 
which enter into most of the models describing the for- 
mation of Schottky barriers [3] . For this reason, theory of 
MIGS were deeply studied in semiconductor heterostruc- 
tures [2, 4-8]. More recently, the attention has turned 
to metal/oxide interfaces, mostly due to several works 
that have first predicted [9, 10] and later demonstrated 
giant tunnel clcctroresistance [11-14] and magnetoresis- 
tance [15-17] in ferroelectric/metal junctions, where the 
decay length of the MIGS plays a major role. 

The eigenstates of the Hamiltonian for the interface 
can be described as the matching at the junction of two 
wave functions: an ordinary bulk Bloch state on the 
metal side, and an exponentially decaying function on 
the insulator side. Assuming that the interface is peri- 



odic in the plane parallel to the boundary [referred to as 
the (x, y) plane], then the components of the wave vec- 
tors parallel to the interface, fey, are real and have to 
be preserved when the electron crosses the junction [18]. 
Therefore, the previous matching is possible only if the 
two wave functions do have the same associated energy, 
symmetry, and fe|| [19]. 

The exponential tails within the insulator decay only 
in the direction perpendicular to the junction (referred 
to as the ^-direction) , and can be actually regarded as 
Bloch functions of the bulk insulator with an associ- 
ated complex wave vector in the perpendicular direction 
k± = k z + iq [18]. The decay length of the wave func- 
tion is given by 1/q. For this reason, despite being a 
genuine interface property, characteristics of MIGS are 
often discussed in terms of the complex band structure 
(CBS) of the bulk insulating material [2]. This analy- 
sis is useful because it allows to predict the character- 
istics of MIGS and other interface properties from the 
simulation of a bulk material, avoiding the much more 
computationally demanding simulation of a realistic in- 
terface. Accordingly, many studies have taken advantage 
of the relatively straightforward interpretation provided 
by the CBS to shed light into interface phenomena such 
as tunneling magnetoresistance [15], or the confinement 
of two-dimensional electron gases [20], with the discus- 
sion often focused on complex bands at high symmetry 
points of the Brillouin zone [21]. 

Despite the previous efforts, a direct comparison of the 
spatial and energetic distribution of the evanescent states 
coming from the CBS (i.e. considering only bulk proper- 
tics) , with the MIGS obtained in a realistic simulation of 
the interface (i.e. taking into account atomic details at 



the junction) are only available in simple molecular elec- 
tronic systems [22]. Up to our knowledge, this study is 
absolutely absent in complex oxide interfaces, where ex- 
tra degrees of freedom such as the spontaneous polariza- 
tion in the insulator might alter and tune their properties. 
We have performed CBS simulations of a bulk insulator 
(PbTiOs) and give a recipe to estimate the decay fac- 
tor of the MIGS from them. The agreement between the 
penetration length estimated from the CBS and the one 
observed after simulations of realistic capacitors is very 
good, especially remarkable when a noble metal (Pt) is 
used as electrode. When the electrode is replaced by a 
metal oxide (SrRu0 3 ), the dependency of the decay fac- 
tor with the the energy shows some extra structure due 
to interfacial effects (symmetry filtering) . In thin film ca- 
pacitors with SrRu03 electrodes, where the polarization 
of the PbTi03 layer is smaller than the bulk spontaneous 
value due to imperfect screening, ferroelectricity does not 
modify strongly the penetration length. Larger values of 
the polarization produce a decrease of the penetration of 
the MIGS. 

We have carried out simulations within the local den- 
sity approximation as implemented in two different codes: 
Siesta [23] for the calculations of realistic interfaces, and 
Quantum-Espresso [24] for the CBS computations of 
bulk materials [25]. The details of the simulations can be 
found in the Supplemental Material [26]. 

Complex band structure of bulk PbTiO^. In order 
to face the difficulty of the problem step by step, we 
have proceed with a centrosymmetric tetragonal phase 
(PA/mmm space group), assuming the in-plane theoret- 
ical lattice constant, a, of an hypothetical^ thick SrTiC>3 
substrate and only allowing the out-of-plane lattice con- 
stant, c, to relax. This is the structure that PbTiC>3 
displays at the center of the unpolarized capacitor. 

In essence a CBS calculation consists in sampling the 
two-dimensional Brillouin zone (2DBZ) parallel to the in- 
terface, and for every fey search for complex fcj_ associated 
with real eigenvalues of the system. Imaginary bands al- 
ways connect extrema of other bands (either real or com- 
plex) [27], and it is often assumed that complex bands 
with the shortest imaginary wave vector q are those con- 
necting the edges of the valence and conduction bands. 
The analysis of the real band structure (see Suppl. Fig. 
2 [26]) reveals that, in the non polar YA/mmm phase, 
PbTi03 possesses a direct gap at X, (%/a, 0, 0), and a 
slightly larger one at Z, (0,0, n/c). 

Complex bands at such high symmetry points have an 
associated k with the form (n/a,0,iq) and (0, 0, ir/c + iq) 
respectively, which correspond to complex bands at fey = 
X and f in the 2DBZ, shown in Fig. 1(a) and (b). We 
observe that the states with the shortest q (largest pen- 
etration lengths) consist of a Ai singlet in the lower en- 
ergy part of the band gap and a A 5 doublet for energies 
close to the conduction band edge. The smaller decay 
factor q takes values from 0.10 to 0.15 in units of (2tt/c) 
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FIG. 1. (color online) Complex band structure of centrosym- 
metric PbTi0 3 at (a) fey = f and (b) fey = X. Black lines 
at central panels correspond to bands with a real value of 
fc_i_. At left and right panels, complex bands with the form 
k± — iq and k± = 7r/c + iq respectively, are plotted as black 
lines. In red we plot complex bands with k z ^ or n/c, with 
the projections on the real plane plotted on the central panel 
and their projections on the imaginary plane on the side pan- 
els (following Ref. [27]). Horizontal dashed lines delimit the 
band gap. In (c) we plot together the imaginary part of bands 
with complex values of k, with fey = f or X. (d) Density of 
states (in arbitrary units) with respect to the energy and the 
imaginary wave vector q. The density of complex bands is 
plotted in a red scale, where darker regions represent a larger 
density of bands. The black solid line represents the effective 
penetration length computed from the CBS (see text). The 
zero of energies corresponds to the top of the valence band. 



at the center of the gap and tends to zero when the en- 
ergy approaches the band edges. In Fig. 1(c) we plot 
together all the energy bands with associated complex 
values of k = (kt\,k z + iq), with fey = T or X. Never- 
theless, we have to keep in mind that the CBS at f and 
X might not be representative of the whole CBS of the 
system, since they are only two points of special sym- 
metry of the total 2DBZ. For a comprehensive analysis, 
we should gather into a plot like Fig. 1(c) the imaginary 
part of complex bands coming from all fey . Instead of 
this, that would result into some of the infinite number 
of complex bands obscuring others, we plot the "density 
of complex bands" in Fig. 1(d). In this 2D-histogram, 
for every value of the energy and imaginary part of the 
wave vector q, we represent in a red scale the number 
of complex bands contributing to that q(E) point. The 



analysis of Fig. 1(d) reveals that the total CBS is clearly 
dominated by bands at q ~* 0.2 and q ~ 0.3. To iden- 
tify the fey that contribute to those clusters of bands we 
resolve the structure and, for every fey , we plot the small- 
est q at a given energy [see Suppl. Fig. 3(a) [26]]. We 
find that from energies between the middle of the gap 
and the bottom of the conduction band, the most pene- 
trating bands come from fen lying in a narrow rectangle 
centered along the f — X path, while from the top of the 
valence band to the center of the gap, the bands with 
smaller q are due to states in small circular (around T) or 
square (around X) regions centered at the high-symmetry 
points. The much larger area of the 2DBZ outside those 
sections contributes to the darker cluster of imaginary 
bands at q ~ 0.2. Similar results have been previously 
reported for other perovskites [15, 21, 28]. The main 
messages that can be drawn might be summarized as 
follows: (i) the most penetrating bands (the most impor- 
tant for tunneling, especially in thick films) come from 
high-symmetry points or lines in the 2DBZ whose rela- 
tive weight is small. That is the reason why those bands 
might be totally inappreciable in the 2D-histogram of 
Fig. 1(d) (see region marked with a square, where the 
bands with q ~ 0.1 at X and T previously discussed were 
expected); and (ii) for ultra-thin dielectric layers or for 
the computation of integrated quantities (such as the in- 
terfacial dipole) it is important to consider higher order 
imaginary bands and non high-symmetry points in the 
2DBZ, since associated phenomena might be dominated 
by them [29]. 

It is sensible to think that the charge transferred from 
the metal to the insulator through the MIGS reflects the 
collective contribution of all complex band in the gap of 
the insulator. Here we also give a recipe to infer from the 
CBS an effective value of the imaginary wave vector for 
This can be estimated by fitting 
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which comparison with the q c ff obtained from the realistic 
simulation of the interface is discussed below. 

MIGS in first-principles simulations of capacitors. The 
explicit influence of the metal/ferroelectric interface in 
the decay of the MIGS is captured in the simulation 
of realistic capacitors from first-principles. Here, we 
have simulated short-circuited (SrRu03) m /(PbTi03) n 
and (Pt)z/(PbTiC>3)p capacitors as model systems, with 
m = 9.5, n = 8.5, I = 9.5 and p = 6.5 unit cells, respec- 
tively. In both cases the PbTiOs films are terminated in 
a Ti02 atomic layer, and we have checked that the band- 
alignment in the most stable paraelectric configuration is 
non-pathological [30]. 

The analysis of MIGS in full capacitor simulations re- 
quires to work with a spatially resolved density of states. 
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FIG. 2. (color online) (a) Schematic representation of the in- 
terface region of the simulated (SrRu03)9.s/(PbTi03)8.5 su- 
perlattice in an unpolarized configuration (Ti atoms in green, 
Ru in magenta, O in blue, Sr in yellow and Pb in gray). Only 
one of the two symmetric interfaces is shown, (b) Layer-by- 
layer PDOS on the atoms at the different layers, as denoted 
by the layer labeling. The plot at the top corresponds to the 
central atomic plane of the PbTiC>3 film. Dashed lines delimit 
one of the energy windows used to perform the integrations 
in Eq. (2). 



Here, we have chosen the layer by layer projected den- 
sity of states (PDOS), as defined in Eq. (21) of Ref. [30]. 
The results for the SrRuG-3 /PbTi03 capacitor in the un- 
polarized configuration arc plotted in Fig. 2(b) (similar 
results are obtained for the Pt/PbTi03 interface). The 
decay of the evanescent tails of the MIGS can be clearly 
observed, with the PDOS showing a clean gap for atomic 
layer further than two unit cells from the interface. 

We can obtain a measure of the spatial distribution 
of the probability density at a given energy integrating 
the PDOS in small windows centered at different energies 
inside the band gap, 



QpDOs{Ei, Zj) 



Ei-- 



P {j,E)dE, 



(2) 



where p(j, E) denotes the PDOS on all the atomic or- 
bitals of a given layer j, located at Zj. The decay of 
this quantity when we move away from the interface pro- 
vides a direct way to obtain the effective decay length 
of the PDOS, 8(E). Since 8 refers to the decay of the 
density, and given the exponential behaviour of the tails 
of the MIGS, we can relate 8 with an effective decay 
of the wave functions as q^ IGS (E) = 1/ [25(E)]. Nu- 
merically, 8(E) is estimated first integrating the PDOS 
after Eq. (2) in energy windows of width AE = 0.02 
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FIG. 3. (color online) (a) Decay of the MIGS for the en- 
ergy window inside the gap of the SrRu03 /PbTi03 capacitor 
shown in Fig. 2(a). The solid curve is a fit to an hyperbolic 
cosine, from which the effective penetration 5(E) is obtained, 
(b) Effective imaginary wave vector q^ s obtained from the 
CBS (red line) and from the fit of the decay of the PDOS in- 
side the insulating material for SrRuOs/PbTiOa (solid black 
line) and Pt/PbTi03 capacitors (dashed black line), (c) To- 
tal DOS of MIGS in the SrRu0 3 /PbTi0 3 (solid black line) 
and Pt/PbTi03 capacitors (dashed black line) obtained in- 
tegrating the layer-by-layer PDOS from the interface to the 
center of the PbTi03 film. In (b) and (c) the zero of energies 
is set to the top of the valence band and the horizontal lines 
indicate the Fermi level of the corresponding capacitor. 

eV. Then, for every energy sampled, we plot the expo- 
nentially decaying Qpdos curve as a function of the dis- 
tance to the interface [Fig. 3(a)], and fit it to a function 
Qpdos(E, z) ~ cosh[z/S(E)]. A hyperbolic cosine was 
used to account for the presence of two interfaces. Re- 
peating this procedure for energy windows covering the 
whole gap we obtain the energy dependence of the effec- 
tive imaginary wave vector, q^- GS (E), that is plotted in 
Fig. 3(b) for the two different capacitors. 

Comparison of decay factors obtained from bulk and in- 
terface simulations. In Fig. 3(b) we plot the effective val- 
ues of the decay factors obtained from a bulk CBS calcu- 
lation [using Eq. (1) and scaling to the band gap obtained 
with Siesta to account for the larger gap obtained with 
Quantum-Espresso], and from the evanescent behavior 
of the MIGS at the interface for PbTi03-based capacitors 
with SrRuOa and Pt electrodes. The agreement between 
q G g S and q^g GS is very good, supporting the idea that 
the penetration of the MIGS is mostly determined by the 
bulk properties of the insulator. The correspondence is 
remarkably good for the capacitor with Pt electrodes, for 
which a featureless dependency of q^ GS with the energy 



is observed. For the metal oxide SrRuOa a more struc- 
tured shape of the decay factor is noticed. In this case, 
to understand the peaks we have to consider interfacial 
effects like the symmetry filtering: the wave functions in 
the electrodes should match evanescent states in the in- 
sulator that are compatible by symmetry. The inspection 
of the bulk band structure of SrRuOa (see Suppl. Fig. 
4 [26]) reveals the presence of a band below the Fermi 
energy with A5 symmetry along the F-X line that nicely 
matches complex bands of PbTi03 of the same symme- 
try, with fc|| along the T — X line of the 2DBZ. Above 
the Fermi energy there are no bands with the appropri- 
ate symmetry to link the complex bands with minimal 
q. The matching has to be done with higher order com- 
plex bands with larger q, explaining the sharp increase 
in <?^ IGS . In the case of Pt, there is a band with the 
appropriate symmetry that crosses completely the band 
gap of PbTi0 3 . 

Finally, not only the decay rate of MIGS is relevant 
but the interfacial DOS of MIGS is also important. We 
find that in this respect the two electrodes behave differ- 
ently as well [see Fig. 3(c)], with Pt/PbTi03 capacitor 
showing a much larger surface density of MIGS than the 
SrRu03/PbTi03 one. This is consistent with the much 
larger DOS in the vicinity of the Fermi level of bulk Pt 
with respect to SrRu03. 

Influence of the ferroelectric polarization. The influ- 
ence of the ferroelectric distortion in the spatial and en- 
ergetic distribution of the evanescent states is known to 
play a mayor role on the electro [21, 28, 31] and magne- 
toresistance [15, 17, 32] of ferroelectric tunnel junctions. 
Unfortunately, a thorough analysis of such effect using 
the arguments developed in this work is limited by band 
alignment issues [30] that restrict the possible electrodes 
and the magnitude of the polarization that can be ex- 
plored. In our case, only the ferroelectric capacitor with 
SrRu03 electrodes and n = 8.5 is non-pathological [30], 
with a polarization within PbTi03 of 24 /iC/cm 2 . Again, 
the agreement between q^ GS and q G g s (computed from 
the CBS of bulk PbTi03 at the same polarization as in 
the capacitor) is very good (see Suppl. Fig. 5 [26]). A 
general increase of the effective decay factor (i.e. shorter 
penetration lengths) with the polarization is observed. 
This growth of q is due to two different factors. First, the 
polarization opens the gap. Second, it pushes up the flat 
band along the T-X line at the bottom of the conduction 
band in the paraelectric configuration (see Suppl. Fig. 
6 [26]). Both facts increase the imaginary part of those 
bands that link to the bottom of the conduction band 
with fen along the T — X line. This effect has already been 
reported for the complex bands at f in BaTi03 [21, 32], 
and was invoked to explain the impact of ferroelectricity 
in the tunnel conductivity of ferroelectric junctions. 

Summary. The main conclusions that can be drawn 
from this work are: (i) although the most penetrating 
bands are located at high-symmetry points, for the study 
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of tunneling across ultra-thin layers or integrated quan- 
tities (where many fey contribute) it is important to con- 
sider complex bands in the whole 2DBZ; (ii) we have 
given a general recipe to estimate from the CBS of the 
bulk insulator the effective decay factor (and its depen- 
dence with the energy within the gap) of the MIGS in 
realistic interfaces, showing a very good agreement with 
the results obtained from first-principles simulations of 
capacitors; (iii) to explain the fine details, especially 
when metal-oxide electrodes are used, a symmetry fil- 
tering analysis has to be performed; and (iv) the ferro- 
electric polarization can be used as a knob to increase 
the decay factor of the CBS. The accuracy of the model 
used to predict the effective decay length of MIGS would 
allow to screen materials to be used as barriers in tunnel 
electro- and magnetoresistance junctions based on knowl- 
edge. We hope this work will encourage experiments and 
set the ground for future theoretical studies of MIGS, 
tunneling and related phenomena. 
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I. SCHEMATIC REPRESENTATION OF METAL-INDUCED GAP STATES (MIGS). 




Suppl. FIG. 1. Schematic representation of metal-induced gaps states in a metal/insulator interface. In (a), the wave function 
of an electronic state with an eigenvalue within the insulator band gap takes the form of a propagating wave function in the 
metallic side that decays exponentially on the insulating side. A schematic representation of the band alignment of the interface 
is shown in (b). MIGS below the Fermi level of the junction are populated and the resulting injection of charge is an important 
contribution to the interfacial dipole that determines the band alignment. 



II. COMPUTATIONAL DETAILS. 



Simulations on this work have been performed within the local density approximation (LDA) in the Ceperley- Alder 
parametrization 1 as implemented in two different codes: Siesta 2,3 for the simulation of realistic capacitors, and the 
Quantum-Espresso ,5 package for the complex band structure calculations. 6 The use of two different codes is a 
delicate issue and thus keeping strict convergence criteria becomes critical in order to ensure the compatibility of all 
the simulations. 

In order to simulate the effect of the mechanical boundary conditions due to the strain imposed by the substrate, 
in all the simulations the in-plane lattice constant was fixed to the theoretical equilibrium lattice constant of bulk 
SrTi0 3 (a — 3.85 A for Quantum-Espresso and a — 3.874 A for Siesta). Since SrRu0 3 is paraelectric at room 
temperature, we have assumed non-spin polarized configurations. 

During the structural optimizations of both bulk PbTiC>3 and the capacitors, a Monkhorst-Pack 7 ' 8 mesh equivalent 
to 12 x 12 x 12 in a perovskite unit cell was used for the sampling of the reciprocal space. Tolerances for the forces 
and stresses were 0.01 eV/A and 0.0001 eV/A 3 , respectively. Other computational parameters, specific to each code, 
are summarized below. 
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A. Quantum-Espresso 

Complex band structure of bulk PbTiC-3 has been calculated using a density-functional plane-wave code with 
pseudopotentials. We have used the ultrasoft Vanderbilt pseudopotentials tabulated in the Vanderbilt Ultra-Soft 
Pseudopotential Site. 9 The plane- wave cutoff was set to 40 Ry. For the generation of the 2D-histogram of Fig. 1(d) 
of the main body of the manuscript, a 48 x 48 mesh of fey points over the two-dimensional Brillouin zone (2DBZ) was 
sampled. 



B. Siesta 

Computations on short-circuited Pt/PbTi03 and SrRu03/PbTi03 capacitors were performed using a numerical 
atomic orbital method, as implemented in the Siesta code. 2 Core electrons were replaced by fully-separable 10 norm- 
conserving pseudopotentials, generated following the recipe given by Troullier and Martins. 11 Further details on the 
pseudopotentials and basis sets can be found in Ref. 12. 

A Fcrmi-Dirac distribution was chosen for the occupation of the one-particle Kohn-Sham electronic eigenstates, 
with a smearing temperature of 8 meV (100 K). The electronic density, Hartree, and exchange-correlation potentials, 
as well as the corresponding matrix elements between the basis orbitals, were computed on a uniform real space grid, 
with an equivalent plane- wave cutoff of 1200 Ry in the representation of the charge density. 

The capacitors were simulated by using a supercell approximation with periodic boundary conditions. 13 A (1 x 1) 
periodicity of the supercell perpendicular to the interface is assumed. This inhibits the appearance of ferroelectric 
domains and/or tiltings and rotations of the O octahedra. A reference ionic configuration was defined by piling up a 
number of unit cells of the perovskite oxide (8.5 for the SrRuC>3-based capacitor and 6.5 for the Pt-based capacitor), 
and 9.5 unit cells of the metal electrode. The PbTiOs films are always terminated in a Ti02 atomic layer. 

To simulate the capacitors in a non polar configuration, we imposed a mirror symmetry plane at the central Ti02 
layer, and relaxed the resulting tetragonal supercells within PA/mmm symmetry. For the ferroelectric capacitors a 
second minimization was carried out, with the constraint of the mirror symmetry plane lifted. We have checked that 
the band-alignment in all systems reported here is non-pathological. 14 



III. BULK BAND STRUCTURE OF PbTi0 3 . 

We compare complex band structure calculations on bulk PbTiC>3 with properties of MIGS in capacitors where the 
out-of-plane cell vectors were allowed to relax. As required by the strong sensitivity of complex bands on the real band 
structure, bulk calculations must be performed under the same symmetry constrains applied to the PbTi0 3 layer in 
the capacitor. It is important to note here that the bulk paraelectric phase discussed in this work is not the cubic 
phase, but a tetragonal centrosymmetric P4/mmm phase. The in-plane lattice constant was fixed to the theoretical 
one of SrTiC"3 and the out-of-plane stress was relaxed while the atomic positions were kept in the centrosymmetric 
positions. 

Given the sensitivity of complex bands to the size of the gap and the curvature of real bands they connect with, 15 
a good agreement between the band structures obtained with both codes is required to obtain comparable results. 
Suppl. Fig. 2 shows that the band structures calculated with the two different codes are virtually identical in shape, 
although the Quantum-Espresso calculation displays a slightly larger gap. This must be taken into account since, 
in principle, a larger gap should translate into slightly larger values of the imaginary part of the complex wave vectors. 



IV. DETAILS ON THE ANALYSIS OF MIGS IN FIRST-PRINCIPLES SIMULATIONS OF 

CAPACITORS. 

The analysis of MIGS in full capacitor simulations requires to work with some sort of energy-resolved probability 
density. For this, a spatially resolved density of states is defined as 

p{i,E)=Y^( dk\(i\ij nk )\ 2 6(E-E nk ), (1) 
where \i) is a normalized function, localized in space around the region of interest. 
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When \i) = |r) is an eigenstate of the position operator, the resulting p(r,E) is commonly known as Zocal density 
of states. However, strong oscillations of this function due to the underlying atomic structure difficult the analysis. 
Alternatively, the nanosmoothed version of this function might be used, but some interfacial properties, like precisely 
the decay length of the MIGS charge in the band gap, are sensitive to the specific convolution function used for the 
nanosmoothing procedure. 16 

A reasonable choice is to work with the layer-by-layer (z- resolved) PDOS (-E-resolved) , defined as in Eq. (1) where 
\i) = \4>nlm) is arL atomic orbital of specified quantum numbers (n,l,m). In this case the bias of the method lies in 
the choice of the basis of atomic orbitals. A sufficiently converged basis should minimize its effect since atoms of the 
same species at different sites are equally described, so the z dependence might be considered as less biased than with 
previous methods. In Fig. 2 of the main body of the manuscript we plot the layer-by-layer PDOS corresponding to 
the PbTi0 3 /SrRu0 3 capacitor. 

The energy distribution of the charge density converges much slower with the fc-point sampling than its spatial 
distribution. For this reason the PDOS was calculated performing an extra non-self consistent calculation with a finer 
fc-point grid of 54 x 54 x 9. 



V. SYMMETRY FILTERING AT THE PbTi0 3 /SrRu0 3 INTERFACE 

Fig. 3(b) of the main body of the text of this work reveals some differences in the properties of MIGS of capacitors 
with SrRu03 or Pt electrodes. The dependence of the effective decay factor with the energy is essentially featureless 
in the case of Pt electrodes, while for SrRuOa this curve presents some structure. The most remarkable characteristic 
in the latter case is a sharp increase in the effective imaginary wave vector in the uppermost part of the PbTi03 gap, 
right above the Fermi level of the capacitor. This behavior is not due to intrinsic characteristics of the complex band 
structure of PbTi03 but to the interface. 

The interface affects the energy dependence of the MIGS penetration mainly through the symmetry filtering: wave 
functions in the electrodes should match evanescent states in the insulator that are compatible by symmetry, filtering 
the complex bands that contribute to the effective value. 

In Suppl. Fig. 3(a) we plot the complex band of minimum q for paraelectric PbTi03 over the 2DBZ. We find 
that from the top of the valence band to the center of the gap, the bands with smaller q are due to states in small 
circular (around T) or square (around X) regions centered at the high-symmetry points, while from energies between 
the middle of the gap and the bottom of the conduction band, the most penetrating bands come from fey lying in a 
narrow rectangle centered along the T — X path. As can be seen in Fig. 1 of the main body of the manuscript, these 
deep-penetrating bands along the high symmetry T — X line are those that link to the bottom of the conduction band 
and have A5 symmetry at the T point. 

The symmetry analysis of the band structure of bulk SrRu03, depicted in Suppl. Fig. 4(a), reveals that there is 
a band (highlighted in red) with the same symmetry as the evanescent states discussed above all along the T — X 
line, and that has its maximum around the Fermi level. For higher energies there are no bands with the appropriate 
symmetry to link the complex bands of minimal q and form deep-penetrating MIGS. The matching has to be done 
with higher order complex bands with larger q. This explains the sharp increase of the effective imaginary part of the 
wave vector for energies above the Fermi level obtained in the PbTiOs/SrRuOs capacitor. 

For Pt there is a band with the appropriate A5 symmetry that crosses completely the band gap, as shown in Suppl. 
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(a) (b) 




Suppl. FIG. 3. (color online) (a) Minimum value of the imaginary part of the complex wave vector (maximum penetration) 
for bulk paraelectric PbTiOs over the 2DBZ. (b) DOS of SrRuC>3 over the 2DBZ. Both quantities are shown for eight different 
values of the energy, referred to the Fermi level of the PbTiOs /SrRu0 3 capacitor. The center of each panel represents the T 
point, while the centers of each side symbolizes the X point. 




Suppl. FIG. 4. (color online) Bulk real band structure of (a) bulk SrRuOs and (b) bulk Pt along the relevant directions in 
the 3-D Brillouin zone of the P4/mmm tetragonal phase. The bands highlighted in red has symmetry A5 along the T-X line, 
and match complex bands in the PbTiOs with the same simmetry with fey along the line F — X. The gray area indicates the 
band gap alignment of the PbTiOs layer for the corresponding capacitors. 



Fig. 4(b). In other words, all the complex bands with minimum q along the T — X path find a symmetry-compatible 
band in the Pt electrode to link with, explaining the featureless shape of q^g GS ■ 

Other interface-intrinsic mechanism that could affect the energy dependence of the q^g is the fey-dependence of 
the DOS of the metal. For instance, for a given energy, a DOS concentrated away from the high symmetry points 
f and X, or the high symmetry path f — X [where the most penetrating complex bands are found, see Suppl. Fig. 
3(a)] would cause an increase of q^g GS at that particular energy. This mechanism does not seem to contribute to the 
sudden increase of q^ 1GS at the Fermi level in the SrRu0 3 /PbTi0 3 capacitor: the DOS of SrRu0 3 varies smoothly 
with the energy and has a significant density around the high symmetry points and paths near the Fermi level [see 
Fig. 3(b)]. 



VI. DISCUSSION ABOUT THE EFFECT OF THE FERROELECTRIC POLARIZATION OF THE 

PbTiOs LAYER 

Most of the discussions in the main text of the manuscript regard the paraelectric phase of the PbTi0 3 /SrRu0 3 
capacitor. However, PbTi0 3 is a ferroelectric perovskite with a large spontaneous polarization in bulk 17 . Indeed, 
the influence of the polarization in the spatial and energetic distribution of evanescent states in ferroelectric/metal 
interfaces has been proposed to contribute largely to the electro- 18-20 and magnetoresistance 21-23 of ferroelectric 
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Suppl. FIG. 5. (color online) (a) Layer-by-layer PDOS of a polar [SrRuOsJg.s/pbTiOajg.s capacitor. The dashed line indicates 
the Fermi level, (b) Profile of the macroscopically-averaged electrostatic potential across the capacitor, (c) Effective imaginary 
wave vectors obtained from the fit of the decay of the PDOS of the PbTi03 layer in the polar capacitor (black line), from 
the CBS of bulk PbTiC>3 with the same ferroelectric distortion as in the center of the capacitor (blue line) and from the CBS 
of centrosymmetric PbTiOs (red line). The two horizontal lines indicate the position of the Fermi level at the two interfaces 
(note that in the capacitor the Fermi level is the same across the whole structure, but to extract q^g GS (E) the PDOS at 
each atomic layer is shifted following the electrictrostatic potential to align the band structures), (d) Complex band structure 
of centrosymmetric P4/mmm (red) and fully relaxed ferroelectric (blue) PbTiOs. The plot was obtained following the same 
recipe as for Fig. 1(d) of the main manuscript. Red an blue solid lines are the effective decays obtained from the CBS of 
centrosymmetric and polarized PbTiOs respectively. Red and blue dashed lines correspond to the bottom of the conduction 
band in each case. The zero of energies in panels (c) and (d) is set to the top of the valence band. 
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Suppl. FIG. 6. (color online) Band structure of bulk ferroelectric PbTiC>3 with the in-plane lattice constant fixed to the 
theoretical one of SrTi03. Band structure calculated with Siesta. 



tunnel junctions. A thorough analysis of such effect, using similar arguments as those developed in this paper, would 
be very useful in this regard. It should be noted, however, that the investigation of those aspects related with the 
realistic interfaces has the limitation imposed by the band alignment issues reported elsewhere 14 , that restrict (i) 
the number of possible electrodes that can be tested within the usual approximations to the exchange-correlation 
functional (only for SrRuC>3, a sizable polarization can be induced in the PbTiOs layer before causing the spurious 
electric breakdown of the capacitor) and (ii) the magnitude of the polarization (even for SrRu03 electrodes, the 
interface becomes pathological for polarizations of the PbTiOs layer larger than ~ 40 /iC/cm 2 ). 

For the capacitors discussed here, we find only the [SrRuOsjg.s/fPbTiOsjs.s capacitor to display a correct band 
alignment in the polar configuration, with a polarization of 24 /iC/cm 2 after relaxing under short-circuit boundary 
conditions. Suppl. Fig. 5(a) displays the layer-by-layer PDOS in the PbTiOs film. The tilting of the bands due to 
the remnant depolarizing field is clearly appreciable. To extract the q^g from the PDOS using the same recipe 
developed in the main manuscript we have to take into account that here the band structure at different atomic 
layers is shifted as a result of the depolarizing field. To perform the fitting of the decay of the MIGS in the same 
fashion as in Fig. 3(a) of the main paper, we align the layer-by-layer PDOS using the macroscopic average of the 
electrostatic potential, shown in Suppl. Fig. 5(b). The resulting q^g GS (E) for the ferroelectric configuration of the 
SrRu03/PbTi03 capacitor is shown as a black line in Suppl. Fig. 5(c). When compared with the q G g S (E) obtained 
for bulk PbTiOs with the same ferroelectric distortion a good agreement is found, as in the paraelectric configuration. 
For such polarization (24 /iC/cm 2 ) the change in the effective decay of MIGS induced by the ferroelectric distortion 
is not very large. 

To appreciate better the effect of the ferroelectricity in the MIGS we compute the CBS for the fully relaxed bulk 
PbTiOs (which displays a polarization of 81 /iC/cm 2 when the in-plane lattice constant is fixed to that of SrTiOs). 
The resulting q™ s is plotted in Suppl. Fig. 5(d). 

The effect of the ferroelectric distortion on the band structure of PbTiOs is two-folded (see Suppl. Fig. 6). In first 
place it opens the gap, tending to increase the imaginary part of the complex bands wave vectors. In addition, it 
pushes up the fiat band that constitutes the bottom of the conduction band along the T-X path in the paraelectric 
configuration. Suppl Fig. 5(d) shows that this increases greatly the imaginary part of those bands that link to the 
bottom of the conduction band with fey along the T — X. Suppl. Fig. 5(d) also shows that, while complex bands 
at high symmetry paths of the 2DBZ arc indeed strongly affected by the polarization of PbTiOs, the dense cluster 
of bands with q ~ 0.2 (2tt/c) display very little distortion in the polar state with respect to the non-polar one. The 
net result is a decrease of the penetration of MIGS, that is larger the closer the energy is to the conduction band, as 
evidenced by the estimated q G ^ s shown in Suppl Fig. 5(d). 
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